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Abstract 

While the Lorentz group serves as the basic language for Ein- 
stein's special theory of relativity, it is turning out to be the basic 
mathematical instrument in optical sciences, particularly in ray op- 
tics and polarization optics. The beam transfer matrix, commonly 
called the ABCD matrix, is shown to be a two- by- two representation 
of the Lorentz group applicable to the three-dimensional space-time 
consisting of two space and one time dimensions. The Jones matrix 
applicable to polarization states turns out to be the two-by-two rep- 
resentations of the Lorentz group applicable to the four-dimensional 
space-time consisting of three space and one time dimensions. The 
four-by-four Mueller matrix applicable to the Stokes parameters as 
well as the Poincare sphere are both shown to be the representations 
of the Lorentz group. 
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1 Introduction 



One complex number contains two independent parameters. One two-by-two 
complex matrix with its four elements contains eight independent parame- 
ters. The unit determinant condition reduces the number of independent 
parameters to six. These matrices form a group which is called the two- 
dimensional special linear group and is denoted as SL{2,c). 

This group has two important subgroups. If we choose only Hermitian 
matrices, this group is SU{2) which is the two-dimensional rotation group 
describing the rotations of electron spins. Its correspondence with the three- 
dimensional rotation group is well known. 

From the group SL{2, c), we can choose matrices with real elements. As 
in the case of the SU (2) subgroup, these two-by-two real matrices contain 
three independent parameters forming the SL{2,r) or Sp{2) subgroups. In 
paraxial ray optics, reflections and refractions are governed by such matrices 
and are called the ABCD matrices (Azzam et al. 1977, Saleh et al. 2007). 
Thus we shall collectively refer these real subgroups as the ABCD group. 

In this report, we first discuss physical applications of this three-parameter 
ABCD group. It is noted that the ABCD matrix can be decomposed into 
three easy-to-understand rotation and squeeze matrices, which can then be 
rotated into a form having equal diagonal elements and two independent 
parameters (Baskal et al. 2009, Baskal et al. 2010). We shall discuss how 
this process can be developed using optical instruments for periodic systems, 
such as laser cavities and multilayer optics. 

While the group SU{2) corresponds to SO{?>) consisting of three-by-three 
matrices, the ABCD group corresponds to the set of three-by-three matrices 
of Lorentz transformations applicable to the three-dimensional space consist- 
ing of one time dimension and two space dimensions. Lorentz transformations 
applicable to z and x directions and rotations around the y axis are governed 
by this group, which is called SO{2, 1). 

Consequently, this aspect of the ABCD group allows us to study the 
essential features of Wigner's little group (Wigner 1939) dictating the internal 
space-time symmetries of relativistic particles moving along the z direction. 
If we add the rotational degree of freedom around the z axis, the group can be 
extended to four-by-four matrices applicable four-dimensional Minkowskian 
space consisting of three space dimensions and one time dimension. 

Using the two-by-two matrix corresponding to the rotation around the z 
axis, we can extend the three-parameter ABCD group to the six-parameter 
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SL{2, c) group and thereafter to the six parameter Lorentz group. In this 
way we can move into polarization optics. 

The physics of hght polarization has a long history (Saleh et al. 2007,Brosseau 
1998). The basic mathematical instruments in this field are the two-by- 
two matrix formahsm for Jones vectors (Jones 1941, 1947) and the four- 
by-four Mueller (Soleillet 1929, Mueller 1943) matrix formalism for the four 
Stokes (Stokes 1852) parameters. For many years, it was known that the 
set of matrices applicable to the two-component Jones vector has the SU (2) 
symmetry. It is shown in this note that this SU{2) symmetry can be ex- 
tended to that of SL{2, c), if we take into account different attenuation rates 
for the two different polarization directions (Opatrny et al.l993, Han et al. 
1997, Ben-Aryeh 2005). 

In this report, we address the issue of unifying these two mathematical 
devices into one mathematical formalism. We show that the Jones vector 
and the Stokes parameters are two-by-two and four-by-four representations 
of the same Lorentz group. Then the question is why we need the Stokes 
parameters with a more complicated four-by-four matrices. The four-by- 
four formalism can tell the degree of coherence between the two orthogonal 
polarization components. However, the Lorentz symmetry cannot change the 
decoherence parameter, in as much as the same way as this symmetry cannot 
change the mass of a given particle. In order to address this issue, we use 
the Poincare sphere. 

The Poincare sphere is a very useful and elegant graphical method to 
represent the polarization state of light (Poincare 1982) In the past it was 
regarded as a geometry of the three-dimensional sphere with a fixed ra- 
dius (Born et al. 1980). This sphere cannot accommodate the symmetry 
of the Lorentz group. Recently, it was noted by the present authors that this 
three-dimensional geometry can be extended to the four-dimensional Lorentz 
group (Ba§kal et al. 2006). It was noted also that this extended Poincare 
sphere, with its Lorentz symmetry, cannot accommodate the variation of 
the decoherence parameter. However, it is possible to introduce two coupled 
Poincare spheres using the S0{?>, 2) deSitter group (Ba§kal et al. 2006). This 
enlarged symmetry allows variations of the decoherence parameter. 

The loss of coherence leads to an increase in the entropy of the system. 
Indeed, the entropy can be computed from this extended Poincare sphere. 
By coupling two Poincare spheres, we can transfer the entropy from one 
sphere to the other. This could serve as another example of Feynman's rest 
of the universe (Feynman 1972). We note here that two coupled harmonic 
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oscillators could serve as an illustrative example of Feynman's rest of the 
universe (Han et al. 1999). 

In Sec. [21 we start with the six generators of the Lorentz group and their 
closed set of commutation relations. Starting from the two-by-two Pauli 
matrices, there are six linearly independent two-by-two matrices. It is noted 
that the three of these matrices are real, and the other three are purely 
imaginary. The ABCD matrix is generated by those imaginary matrices. 

In Sec. [3l it is shown that the optical ABCD matrix can be decomposed 
into a product of three convenient matrices which are rotation and squeeze 
matrices. In Sec. HJ we use these properties to deal with periodic systems in 
optics. 

In Sec. El the mathematical language of the ABCD matrix is translated 
into that of the internal space-time symmetries of relativistic elementary 
particles. We note that the Wigner decomposition and the Bargmann de- 
composition can be explained in terms of the decomposition properties of the 
ABCD matrix. 

In Sec. [HI it is shown that the two-by-two Jones vector formalism is a 
representation of the Lorentz group. We add a squeeze matrix to the con- 
ventional SU{2) formalism. The symmetry group is SL{2,c) corresponding 
to the full Lorentz group with six independent parameters. From this, it 
is possible to construct a four-by-four representation of the Lorentz group. 
In Sec. [71 we show that the Mueller matrix applicable to the four Stokes 
parameters is the four-by-four representation of the Lorentz group. 

Unlike the case of the two-by-two representation applicable to the Jones 
vectors, the four-by-four representation can deal with the decoherence prob- 
lems. On the other hand, the decoherence parameter remains as a Lorentz- 
invariant parameter. In order to deal with this problem, we enlarge the 
symmetry group from the traditional Lorentz group of 5*0(3, 1) to the deSit- 
ter group of 0(3.2). We carry out this operation by extending the concept 
of the Poincare sphere. 

In Sec. [HI we first extend the traditional three-dimensional sphere to a 
four-dimensional object, in order to take into account all the symmetry con- 
tents of the Stokes parameters. We then introduce two Poincars spheres 
coupled to each other through the symmetry of the 0(3,2) deSitter group. 
This allows the exchange of the decoherence parameters between the two 
spheres. Furthermore, this picture of the two coupled Poincare spheres con- 
stitute another example of Feynman's rest of the universe (Feynman 1972). 
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2 Group of ABCD Matrices 



The Lorentz group is the group of four-by-four matrices apphcable to the 
four-dimensional Minkowskian space of {x,y,z,t). The group is generated 
by three rotation generators Jj and three boost generators Ki, satisfying a 
closed set of commutation relations 

[Jj, Jj] — ieijkJk, [Ji: Kj] — ieijkKk, [Ki, Kj] — —ieijkJk- (1) 

These generators are all four-by-four matrices, and we shall give their explicit 
expressions in later sections. 

Fortunately, there is a set of two-by-two matrices which satisfy the same 
commutation relations. We can write them as 

They are the Pauli spin matrices which are very familiar to us. They are 

Hermitian matrices. 

The boost generators can take the form oi iJi, or explicitly 

Unlike Jj, these matrices are anti-Hcrmitian. 

The group of two- by-two matrices generated by the above six two-by-two 
matrices is called SL{2, c). Since the generators of this group share the same 
set of commutation relations as the Lorentz group, they are said to be locally 
isomorphic to each other. In this paper, we shall avoid this mathematical 
word, and simply say that SL{2, c) is "like" the Lorentz group. 

There are a number of interesting subgroup of this SL(2, c) group. The 
three generators Jj generates the SU{2) subgroup which is like the three- 
dimensional rotation group. This aspect is well known. Three Ki alone do 
not form a closed set of commutation relations. 

Among the six generators, J2,Ki, and K^ are pure imaginary, and they 
generate two-by-two matrices with real elements. Furthermore they satisfy 
the closed set of commutation relations 

[J2, Ki] = -iKs, [J2, Ks] = iK,, [Ki, Ks] = iJ^. (4) 

The group of two-by-two matrices generated by these three matrices is called 
Sp{2) or the two-dimensional symplectic group. It is like the Lorentz group 
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applicable to the three-dimensional space of {z,x,t). However, this group 
constitutes the fundamental mathematical language for the optical beam 
transfer matrix often called the ABCD matrix. 

The ABCD matrix is a two- by-two matrix with real elements, and its 
determinant is one. There are therefore three independent parameters. These 
elements are determined by optical materials and how they are arranged. 
The purpose of this note is to explore its mathematical properties which can 
address more fundamental issues in physics. 

First of all, the trace of this matrix could be less than two, equal to two, 
or greater than two. We are interested in what physical conclusions we can 
derive from these numbers. 

In order to bring the ABCD matrix to the form which will describe 
the particle symmetries, we should first transform it into the equi-diagonal 
form where the two diagonal elements are equal to each other (Ba§kal et 
al. 2009, Ba§kal et al. 2010). We can achieve this goal by a similarity 
transformation with a rotation matrix. Thus, this rotation angle becomes one 
of the three independent parameters, and the equi-diagonal ABCD matrix 
has two independent parameters. 

We shall call this equi-diagonal matrix the core of the ABCD matrix, 
and use the notation [ABCD]. This matrix is not always diagonalized. This 
creates non-trivial problem. We shall examine how optical devices, especially 
periodic systems, can lead us to a better understanding of the problem. For 
this purpose, we discuss laser cavities and multilayer systems in detail. 

If the trace is less than two, the core can be written as 

\ABCn]-( -e''sin(7/2)\ 
^^^^^J- U-''sin(7/2) cos(7/2) J' 

The diagonal elements are equal and smaller than one. 
If the trace is greater than two, the core takes the form 

\ARCr,\-( cosh(7/2) e''sinh(7/2)\ 
^ ^~Ve-''sinh(7/2) cosh (7/2 ) ' 

Here again the diagonal elements are equal, but they are greater than one. 
If the trace is equal to two, the [ABCD] matrix becomes 

[ABCD]^{^1 7). (7) 

This matrix also has the same diagonal element, and they are equal to one. 
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The triangular matrix of Eq. ([7]) cannot be diagonalized. As for the ma- 
trices of Eq.([S]) and Eq.([n]), they can be diagonahzed, but not by rotation 
alone. These mathematical subtleties are not well known. The purpose of 
this report is to show how much physics we can understand by studying this 
mathematical subtleties. 

The mathematics of group theory allows us to write down a four-by-four 
Lorentz-transformation matrix for every two-by-two matrix discussed in this 
paper. In this way, the three matrices given in Eq.([5]), Eq.([6]), and Eq.(I7]) 
lead to the internal space-time symmetries of elementary particles. They 
respectively correspond to the symmetries of massive, imaginary-mass, and 
massless particles, respectively (Wigner 1939, Kim et al. 1986). 



3 Decomposition of the ABCD Matrix 

We are interested in writing the three different forms of the core matrix in 
one expression. 

where the parameters x and y are determined by the optical materials and 
how they are arranged. The exponent of this matrix is 



1 / —X — y 
2\x-y 

li X > y, the exponent becomes 

7/0 — exp (77) 
2\exp{-ri) 

which leads to the core matrix of Eq. ([5]) with 



(9) 



(10) 



7 = - 2/^ 



e" = (11) 

\ x-y 

The core matrix [ABCD] can be written as a similarity transformation 

[ABCD] = B{7])R{e)B{-7]) (12) 
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with 

Biv) 

Rid) 



ev/2 



e-''/^ 



cos(0/2) -sin(^^/2)\ 

sm{e/2) cos(^/2) r ^ ^ 



where 7 is now replaced by the rotation angle 6. R{0) is a rotation matrix, 
and B{ri) is a squeeze matrix. 
If a; < the exponent becomes 



7 / — exp (77) 

2V-exp(-?7) 

leading to the core matrix of Eq. ([6]) , with 



(14) 



7 = \/y^ -x^, 



= f-^. (15) 

The matrix can now be decomposed into a similarity transformation 

[ABCD] = B{r])S{-X)B{-r]), (16) 

with 

S(A)=("*(V2) .„h(A/2U 

V sm(A/2) cosh(A/2) y ^ ^ 

where 7 is replaced by the boost parameter A. The matrix B{ri) takes the 
diagonal form given in Eq. (fT2|) with 77 defined in Eq. (fT5|l . S{\) is a squeeze 
matrix. 

If X = y, the exponent becomes 



1 —X 

1 



(18) 



with X = y = •y. 

We now have combined three different expressions for the core of the 
ABCD matrix into one exponential form of Eq.([8]). This form can be de- 
composed into three matrices constituting a similarity transformation. 
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There is another form of decomposition known as the Bargmann decom- 
position (Bargmann 1947), which states that the core of the ABCD matrix 
can be written as 

[ABCD] = R{a)S{-2x)Ria), (19) 

where the forms of the rotation matrix R and the squeeze matrix S are 
given as in Eq.f lT^ and Eq.f lT7|) respectively. If we carry out the matrix 
multiphcation, the [ABCD] matrix becomes 



(cosh x) cos a — sinh x — (cosh x) sin a 

— sinh X + (cosh x) sin a (cosh x) cos a 



(20) 



This matrix also has two independent parameters a and x- We can write 
these parameters in terms of 6 and r] by comparing the matrix elements. For 
instance, if x > y, the diagonal elements lead to 

cos(^^/2) = (coshx) cosa. (21) 

The off-diagonal elements lead to 



2^ (cosh x) sin a + sinh x 
(cosh x) sin a — sinh x 



(22) 



As for physical applications, let us consider periodic systems, such as laser 
cavities and multilayer systems. The exponential form given in Eq.® tells 
us that it is a matter of replacing the 6 parameter by NO for N repeated 
applications (Ba§kal et al. 2010). Let us see some examples. 



4 Periodic Systems in Ray optics 

Let us summarize the content of Sec. [31 First of all, the optical ABCD matrix 
can be brought to an equi-diagonal form by a similarity transformation, and 
we call this equi-diagonal matrix the core of ABCD, and use the notation 
[ABCD]. Thus, 

(ABCD) = T [ABCD] T-\ (23) 

where we use {ABCD) for the original ABCD matrix. The transformation 
matrix can be a rotation or a triangular matrix. In addition, the core matrix 
can be written as a Wigner decomposition of the form 

[ABCD] = B{r]) 1^(7) B{-r]), (24) 



9 



which is another similarity transformation with B{ri) given in Eq.f ll3p . The 
W{'y) is one of the three matrices 

/cos(7/2) -sin(7/2)\ / cosh(7/2) -sinh(7/2)\ /I -7\ 
Uin(7/2) cosh(7/2)j' ^ - sinh(7/2) cosh(7/2) j ' [o ij' 

(25) 

and 

(ABCD) = [T Bir^)] Wi^) [T Biv)]-\ (26) 

For repeated apphcations of the ABCD matrix we need an expression 
of the form (ABCD)^ . Thanks to this form of similarity transformation, 
(ABCD)'^ is now simplified as 

(ABCD)^ = [T B{ri)] W{N-f) [T B{ri)]-\ (27) 

Thus the problem of periodic systems in optics is to find the core matrix 
[ABCD]. 



4.1 Laser Cavities 

As the first example, let us consider the laser cavity consisting of two identical 
concave mirrors separated by a distance d. Then the ABCD matrix for a 
round trip of one beam is 



1 0\/l d\/ 1 0\ /I d 
-2/R l) \0 l) \-2/R l) \0 1 



where the matrices 



1 0\ fl d 
-2/R ij' [o 1 



(28) 



(29) 



are the mirror and translation matrices respectively. The parameters R and 
d are the radius of the mirror and the mirror separation respectively. This 
form is quite familiar to us from the laser literature (Yariv 1975, Haus 1984, 
Hawkes et al. 1995). 

However, the main issue here is how to alleviate the problem of taking 
the Nth power of chains of matrices which corresponds to the repetition of 
this process for many times. This can be achieved when the matrices in 
Eq. fl28l) can be brought to an equi-diagonal form and eventually to a form 
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of the Wigner decomposition. Thus, we are interested in finding the core of 
Eq.f l28p . For his purpose, we rewrite the matrix of Eq. fl28p as 



1 d/2\ f 1 
1 ) \-2/R 1 



(30) 



1 OWl 
-2/R Ij \0 

In this way, we translate the system by ~d/2 using a translation matrix 
given in Eq.f l29p . and write the ABCD matrix of Eq. (!28|) as 



-d/2 
1 



l-d/R d-d^/2R\ 



d/2 
1 



-2/R l~d/R ) 

We are thus led to concentrate on the matrix in the middle 

l-d/R d-d^/2R\ 
-2/R l-d/R y ' 

which can be written as 

fVd \fl-d/R l-d/2R\/l/Vd 
[o l/Vd)\-2d/R l-d/R J [ Vd 

It is then possible to decompose the ABCD matrix into 

E E-\ 



(31) 



(32) 



(33) 



(34) 



with 



C 



E 



l-d/R l-d/2R\ 
-2d/R l-d/R j ' 

1 -d/2\ /Vd 
1 j V 1/Vd 



(35) 



The C matrix now contains only dimensionless numbers, and it can be written 
as 

cos(7/2) e''sin(7/2)\ 
-e-^ sin(7/2) cos(7/2) ) ' 

with 



C 



(36) 



cos(7/2) 



1-^. 

R' 



l2R-d 
Ad ■ 



(37) 
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Here both d and R are positive, and the restriction on them is that d be 
smaller than 2R. This is the stability condition frequently mentioned in the 
literature (Haus 1984, Hawkes et al. 1995). 

Thus, the [ABCD] core matrix is C^, and takes the form 

\Afirn]-f e^sin(7)\ , . 

^^^^^J- Ue-''sin(7) cos(7) j' ^^^^ 

and the similarity transformation which connects this core matrix with the 
original ABCD matrix of Eg. (1281) is E given in Eg. (1351) . 



4.2 Multilayer Optics 

We consider an optical beam going through a periodic medium with two 
different refractive indices. If the beam traveling in the first medium hits the 
second medium, it is partially transmitted and partially reflected. In order 
to maintain the continuity of the Poynting vector, we define the electric fields 

as 

Ei^^ = exp (±ikiz — ut) 

E^f^ = ^= exp {±ik2Z - cut) (39) 



n2 

for the optical beams in the first and second media respectively. The super- 
script (+) and (— ) are for the incoming and reflected rays respectively. 

These two optical rays are related by the two- by-two ABCD matrix, 
according to 

Of course the elements of this matrix are determined by transmission coeffi- 
cients as well as the phase shifts the beams experience while going through 
the media (Azzam et al. 1977, Georgieva et al. 2001). 

When the beam goes through the first medium to the second, we may use 
the the boundary matrix given in (Azzam et al. 1977) and in (Monzon et al. 
2000, Monzon et al.2002). In terms of the refractive indexes rii and n2, this 
matrix is 

- sinh((T/2)\ , . 

Uinh((T/2) cosh((T/2) j ' ^ ^ 
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where one can write the a parameter as 



cosh [^-] = ?i±^, sinh C"-] = Ji^. (42) 



2/ 2^Jn^n^ \2 J 2yjn^n^ 

The boundary matrix for the beam going from the second medium should be 

Q{-a). 

In addition, we have to consider the phase shifts the beams have to go 
through. When the beam goes trough the first media, we can use the phase- 
shift matrix 

^(^i)=( e^'/O' ^^^^ 

and a similar expression for P {82) for the second medium. The phase shift 
5 is determined by the wave number and the thickness of the medium. 

We are thus led to consider one complete cycle starting from the midpoint 
of the second medium, and write 

P {82/2) Q{a)P (5i) Q{-a)P {82/2) . (44) 

There are two questions in regards to the above matrix multiplication. 
One is whether each matrix in this expression can be converted into a matrix 
with real elements and the other is whether the result of this matrix multi- 
plication accommodates Wigner and Bargmann decompositions. In order to 
answer the first question, let us consider the similarity transformation 

P{8)Q{a) (45) 

with 



This transformation leads to 



R{8)Q{a), (47) 



where 

•cos((5/2) -sin((5/2)\ 



This notation is consistent with the rotation matrices used in Sec. [31 
Let us make another similarity transformation with 

--71 (A 
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This changes Q{(j) into B{a) without changing R{S), where 

again consistent with the B{ri) matrix used in Sec. [31 

Thus the net similarity transformation matrix is (Georgieva et al. 2001) 



1 / e*'^/^ e*'^/'^ \ 
C = C2C1 = _g-W4 ) • (51) 



If we apply this similarity transformation to the long matrix chain of 
Eq. (H4l) . it becomes another chain 

M = R (52/2) Bia)R ((5i) 5(-a)i? (52/2) , (52) 

where all the matrices are real. 

Let us now address the main question of whether this matrix chain can 
be brought to one equi-diagonal matrix. We note first that the three middle 
matrices can be written in a familiar form 

M = B{a)R{6i)B{-a) 

cos(5i/2) -e'^ sin(5i/2) \ 
e--^ sin(5i/2) cos(5i/2) J ^ ^ 

However, due to the rotation matrix R{62/2) at the beginning and at the 
end of Eq. (!52!) . it is not clear whether the entire chain can be written as a 
similarity transformation. 

In order to resolve this issue, let us write Eq. 0531) as a Bargmann decom- 
position 

R{a)S{-2x)Ria), (54) 

with its explicit expression given in Eq. (l20l) . The parameters a and x ^''^^ 
related to a and 61 by 

cos(5i/2) = (coshx) cosa, 

9^ (cosh y) sin a + sinh y . . 

^ = / ■ -f-^- 55 

(cosh X) sm a — smh x 
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It is now clear that the entire chain of Eq. f l44p can be written as another 
Bargmann decomposition 

M = R{a + 52/2)S{-2x)R{a + 52/2). (56) 

Finally, this expression can be converted to a Wigner decomposition (Georgieva 
et al. 2003) 

M = B{r])R{e)B{-r]), (57) 

with 

cos(^/2) = (coshx) cos(a + 52/2), 

^2ri ^ (cosh x) sin(a + 62/2) + sinh x /^gx 
(coshx) sin(a + ^2/2) — sinhx 

The decomposition of Eg. (1571) allows us to deal with the periodic system of 
multilayers. For repeated application of M, we can now write 

= B{r])R{Ne)B{-7]). (59) 



5 Space-time Symmetries 

In Sec. [3l we have seen that the two-by-two matrices can provide a very pow- 
erful language for optical systems. However, this language is not restricted 
to the two-dimensional world. It can be translated into the four-dimensional 
world of Einstein's special relativity where Lorentz transformations play the 
central role. 

In mathematics, the group of two-by-two unimodular matrices is called 
SL{2,c). The group of four-dimensional matrices performing Lorentz trans- 
formations on the Minkowskian four-vector {t,z,x,y) is called the 50(3,1) 
Lorentz group. The group SL{2,c) has six generators, so does the 50(3, 1), 
but the corresponding matrices are two-by-two and four-by-four, respectively. 
Their generators satisfy the same set of commutation relations as in Eq.([T]). 
This correspondence is called the local isomorphism between the SL{2, c) 
and 50(3, 1) groups. The four- dimensional generators are 



Ji = 



/o 








0^ 



































i 


, J2 = 








—i 


































i 

















—i 


\o 


—i 





0) 




[0 








0) 




\o 





i 


; 



(60) 
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and 



/O z ON^ 



1 
VO 0/ 





1 

Vo oy 



/O z ON^ 

z 



Vo oy 



(61) 

This mathematical property allows us to explain events in Einstein's 
Lorentz-covariant world in terms of what we observe in optics laboratories. 



5.1 Two-by-two and four-by-four representations of the 
Lorentz group 



The content of this correspondence is somewhat complicated, but for the 
present purpose, we can start with the Minkowskian four-vector (t, z, x, y) 
written as 

t + z X — iy' 
+ iy t — z 



X 



(62) 



whose determinant 



j2 2 2 2 

t — z — X — y 



(63) 

is left invariant under Lorentz transformations. Now consider a transforma- 
tion 

X' = G X G\ (64) 

where G a unimodular matrix whose determinant is one. Let us write this 
matrix as 

'a P' 
,7 ^ . 

where the elements can be complex numbers. If the determinant of this 
matrix is one, there are only six independent parameters. Thus this matrix 
can be generated by the six generators given in Eq.([2]) and Eq.([3]). 
The transformation of Eq.f l64p can be explicitly written as 



G 



(65) 



t' + z' 
x' + iy' 



X — ly 
t'-z' 



t + z 
X + iy 



X ■ 
t 



ly 

■ z 



(3* 



6* 



(66) 



We can now translate this formula into 



f t' + z' \ 
x' — iy' 
x' + iy' 

\ t'-z' J 



/ aa* a/3* Pa* 

0:7* a6* PY 

7a* 7/3* 6a* 

V 77* 75* 57* 



PP*\ / t + z \ 
P6* X - iy 
6P* X + iy 
65* J \t-z J 



(67) 
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It can be seen that the above transformation matrix can be expressed in the 
form of a Kronecker product as 



G®G* 



where 



G* 



(3* 
5* 



(68) 
(69) 



Then the components of X' are related to (t', z', x', y') as 







/I 


1 








\ 


f t' + z' \ 


z' 


1 


1 


-1 










t' -z' 


x' 


~ 2 








1 


1 




x' — iy' 


w) 




\o 





i 


—i 


J 


\x' + iy' / 



(70) 



Likewise, the two-by-two matrix for the four-momentum of the particle 
takes the form 

PO +Pz Px- iPy ' 
.Px + iPy Po-Pz 



P 



(71) 



with po = ^m? + pI + Px+ p1- The transformation of this matrix takes the 
same form as that of the space-time four- vector given in Eq. (16^ and Eq. (!66|) . 



5.2 Internal Space-time Symmetries of Elementary Par- 
ticles 

These properties are applicable to many other branches of physics. For in- 
stance, one of the persisting problem is the internal space-time symmetry of 
elementary particles in Einstein's Lorentz-covariant world. The mathematics 
of group theory allows us to translate the rotation and squeeze matrices of 
Eq.( |T3i) and Eq.(fT71) into the following four-by-four matrices respectively. 



R{e) 



/I 













COS 9 


— sin 9 








sin 6* 


cos 9 





Vo 








1/ 



5(A) 



/ cosh A sinh A \ 

10 

sinh A cosh A 

V 1/ 



( cosh Tj 

sinh 7] 


V 



sinh 7] 
cosh 7] 









1 




0\ 





1/ 



(72) 
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z{4>) 



(73) 



They are applicable to the Minkowskian four- vector {x,y,z,t). The R{0) 
matrix performs a rotation around the y axis, and S{X) is for Lorentz boosts 
along the x axis. The B{ri) matrix boosts the system along the z direction. 

Together with a rotation matrix around z axis (Han et al. 1986) 

/I \ 

10 

COS (p — sin 4 

\0 sin0 COS0 / 

they constitute Wigner's little groups dictating internal space-time symme- 
tries of massive and imaginary-mass particles (Wigner 1939). The triangular 
matrix of Eq.(I71) leads to the little group for massless particles. The little 
groups are the subgroups of the Lorentz group whose transformations leave 
the four- momentum of a relativistic particle invariant. 

It is possible to compute the two-by-two equivalent of the above Z{(f)) 
matrix using the relation given in Eq.f l64p . It takes the form 

This matrix contains complex elements. This is the reason why it is not 
mentioned in our discussions of the ABCD matrix. This rotation matrix 
will play an important role in polarization optics which will be discussed in 
Sec. Eland Sec. [71 

Let us go back to Eq.(|5]) which, according to Eq.(IT2l). can be decomposed 
to a similarity transformation 

W{t], 9) = B{r])R{e)B{-7]). (75) 

We can write this decomposition with the four-by-four matrices given in 
Eq.(I72D. 

Let us then consider a massive particle moving along the z direction with 
the velocity parameter v/c = tanh?7, and its four-momentum 

(m cosh?], m sinh?], 0, 0), (76) 

where m is the mass of the particle. 

We can boost this particle using the boost matrix B{—ri), which is the 
inverse of the four-by-four matrix given in Eq.f l72|) . The particle becomes at 
rest, with its four- momentum 

(m, 0,0,0), (77) 
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Figure 1: Wigner decomposition (left) and Bargmann decomposition (right). 
These figures illustrate momentum preserving transformations. In the 
Wigner transformation, a massive particle is brought to its rest frame. It 
can be rotated while the momentum remains the same. This particle is 
then boosted back to the frame with gaining original momentum. In the 
Bargmann decomposition, the momentum is rotated, boosted, and rotated 
to its original state. 



and with zero velocity. The rotation matrix R{0) rotates this particle with- 
out changing its momentum. During this process, the particle changes the 
direction of its spin. Finally, B{r]) boosts the particle and restores its mo- 
mentum. In this way, the four- by-four expression for Eq. (fT2l) changes the 
internal space-time structure of the particle. 

One key question from this table is what happens to the 0(3)-like little 
group when the particle momentum becomes infinity or its mass becomes 
zero. The question is whether the little group for a massive particles becomes 
that of a massless particle. The answer to this question is Yes, but this issue 
had a stormy history before this definitive answer (Kim et al. 1990). Indeed, 
when 7] becomes infinity, the four-by-four form of Eq flTSl) becomes 



/1 + 7V2 


7V2 
1 + 7V2 


7 




7V2 


7 





7 


7 


1 





V 7 


7 





1/ 



(78) 



When applied to the momentum of a massless particle moving in the negative 
z direction with 

(p,-p,0,0), (79) 
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it leaves the above four- momentum invariant, but on the other hand performs 
a gauge transformation when apphed to the four-potentials of the electromag- 
netic field (Kim et al. 1990). 



6 Jones Vectors 

In studying polarized light propagating along the z direction, the traditional 
approach is to consider the x and y components of the electric fields. Their 
amplitude ratio and the phase difference determine the state of polarization. 
Thus, we can change the polarization either by adjusting the amplitudes, 
by changing the relative phases, or both. For convenience, we call the opti- 
cal device which changes amplitudes an "attenuator" and the device which 
changes the relative phase a "phase shifter." 

The traditional language for this two-component light is the Jones-matrix 
formalism which is discussed in standard optics textbooks (Hecht 1970). In 
this formalism, the above two components are combined into one column 
matrix with the exponential form for the sinusoidal function 



/ y4exp {2(^2; — Ctjt -|- </)i)} 
Ey) V-B exp {i(A;2; — -|- 02)} 



(80) 



This column matrix is called the Jones vector (Jones 1941) . 

The Jones-matrix formalism starts with the projection operator (Hecht 
1970) 

\ :). (-) 

applicable to the Jones vector of Eq. flHOj) . This operator keeps the x compo- 
nent and completely eliminates the y-component of the electric field. 

This is an oversimplification of the real world where the attenuation factor 
in the y direction is greater than that of the x direction. We shall replace 
this projection operator by an attenuation matrix which is closer to the real 
world. 

In this note, we replace the projection operator of Eq.( l8T]) by a squeeze 
matrix. There are two transverse directions which are perpendicular to each 
other. The absorption coefficient in one transverse direction could be different 
from the coefficient along the other direction. Thus, there is the "polariza- 
tion" coordinate in which the absorption can be described by (Opatrny et 
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al. 1993, Han et al. 1997, Ben-Aryeh 2005) 



with /i = /i2 — /xi . Let us look at the projection operator of Eq. llHTl) . Physi- 
cally, it means that the absorption coefficient along the y direction is much 
larger than that of the x direction. The absorption matrix in Eq.( l82p be- 
comes the projection matrix if /ii is very close to zero and /i2 becomes in- 
finitely large. The projection operator of Eq. flHT]) is therefore a special case 
of the above attenuation matrix. 

The attenuation matrix of Eq. fl5^ tells us that the electric fields are at- 
tenuated at two different rates. The exponential factor e"^'^^"''^^'*/^ reduces 
both components at the same rate and does not affect the state of polariza- 
tion. The effect of polarization is solely determined by the squeeze matrix 

/ e''/^ \ 
^(/^)=(o e-'/O' ^^^^ 

which is given in Eq.f lT3|) . This type of mathematical operation is quite 
familiar from studies of squeezed states of light, if not from Lorentz boosts 
of spinors. Thus, we are expanding the Jones-matrix formalism by replacing 
the projection operator of Eq. flST]) by the squeeze operator in Eq.f lH^ . 

Another basic element is the optical filter with two different values of the 
index of refraction along the two orthogonal directions. The effect on this 
filter can be written as 



e''^ J'"" I e''/^ ^ ' 

with 5 = 81—82 ■ In measurement processes, the overall phase factor e*(''i+'^2)/2 
cannot be detected, and can therefore be deleted. The polarization effect of 
the filter is solely determined by the matrix 

^W=( e-'/O' ^^^^ 

which leads to a phase difference of 8 between the x and y components. The 
mathematical expression for this matrix is given in Eq. ( FM|) . It has a different 
physical meaning in the symmetry of the Lorentz group. 
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The polarization axes are not always the x and y axes. For this reason, 
we need the rotation matrix 



Vsin(t^/2) cos(t^/2) / ^ ' 

The traditional Jones-matrix formalism consists of systematic combinations 
of the above three components given in Eq. flST]) . Eq. flH^ and Eq.( 



6.1 Squeeze and Phase shift 

The effect of the phase shift matrix Z{5) of Eq. (!85!) on the Jones vector is 
well known, but the effect of the squeeze matrix of Eq. (l5^ is not addressed 
adequately in the literature. Let us discuss the combined effect of these two 
matrices. First of all both are diagonal and they commute with each other. 

The effect of the squeeze matrix on the Jones vector is straight-forward 
and the net result is 



e-^'l^ )\Ey) ~ W^'^Ey 



(87) 



This squeeze transformation expands one amplitude, while contracting the 
other so that the product of the amplitudes remain invariant. This squeeze 
transformation is illustrated in Fig. [2j 

In order to illustrate phase shifts, we start with the Jones vector of the 
form 

exp {ikz) 
,exp [i{kz — vr/2)] 



whose real part is 

x\ _ f cos(fcz) 
,y) Vsin(fcz) 

which corresponds to a circular polarization with 



(88) 
(89) 



x^ + y'^ = l. (90) 
If we apply the phase shift matrix, the resulting vector is 

x \ _ ( cos(fcz + 5/2)\ , , 

y) \sin{kz-5/2))' ' 
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Phase Shift 



Squeeze 




Figure 2: Squeeze and Phase Shift. Both squeeze and phase shifts resuh in 
eUiptic deformations, but they are done differently. 



which can be written as 



x\ _ f cos{kz — tt/A + a) 
yj \cos{kz — tt/A — a) 



with 



Then 



and 



(92) 



« = ^ + ^- (93) 



X + y = 2(cosa) cos{kz — 7r/4), 

a; — ?/ = — 2(sin a) sin(A;z — 7r/4), (94) 



{x + y)^ ^(x-yY (95) 



4(cosa)^ 4(sina)^ 

This is an elhptic polarization. 

The squeeze operation of Eq. flSSjl is relatively simple. It changes the am- 
plitudes, and it commutes with the phase shift matrix. Thus, the combined 
effect could be illustrated in Fig. |2l 

6.2 Rotation of the Polarization Axes 

If the polarization coordinate is the same as the xy coordinate where the 
electric field components take the form of Eq. flHOj) . the attenuator is directly 
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applicable to the Jones vector as in Eq. fl87p . If the polarization coordinate 
is rotated by an angle of (6/2), or by the matrix 

then the polarization coordinate makes an angle (0/2) with the xy coordinate 
system. So the phase shifter takes the form 



Z{e,S) = Ri9)P{S)R{-9) (97) 

(98) 



cos((5/2) + i sin(5/2) cos 6 i sin(5/2) sin 6 

i sin(5/2) sin6' cos(5/2) — i sin(5/2) cos 6* 

If the polarization coordinate system is rotated by 45°, the phase shifter 
matrix becomes 

In order to illustrate what this matrix does to the polarized beams, let 
us start with the circularly polarized wave 



1 

—i 



^{ikz-iu^t)^ (100) 



whose real part is 

X\ _ f cos(kz — ujt) 
Y J \ sm(kz — cot) 

This leads to the familiar equation for the circle 



(101) 



X^ + Y^ = l. (102) 



If the phase shifter of Eq.( !99l) is applied to the above Jones vector, the 
result IS 

/ [cos(5/2) + sin(5/2)] cos{kz - cot) \ , . 

\i[sm{S/2) - cos(5/2)] sm{kz - cot) ) ^ > 

with 

cos(5/2) = cos {[5/2 + 7r/4] - 7r/4) , 

sin(5/2) = cos {[5/2 + 7r/4] + tx/A) . (104) 

Thus, 
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cos(5/2) + sm(5/2) = ^2 cos (5/2 + 7r/4) , 
cos(5/2) - sin(5/2) = ^2 sin {5/2 + 7r/4) . (105) 
After the phase shift, the Jones vector becomes 



/ [a/2 cos a] cos{kz — cot) 
V [v^sina] sin(A;2; — cot) 



(106) 



with 

« = ^ + ^- (107) 
The the x and y components will satisfy the equation 

+ = 1- (108) 

(y2cosa)2 (y2sina)2 ^ ^ 

This is an elliptic polarization. 

Let us next consider rotations of the squeeze matrix 

B{e, 12) = R{e)B{i2)R{-e), (io9) 

which leads to 

'cosh(/i/2) + sinh(/i/2) cos^ sinh(yu/2) sin6' ^ 

) cos^ 

(110) 



^^^'^^^ ^ sinh(/i/2) sin e cosh(/i/2) - sinh(^/2) cos^. 
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From Sec. |31 we are familiar with this squeeze operation which changes the 
amphtudes. 

If two squeeze transformations are made in two different directions, the 
result is not another squeeze, but a squeeze matrix followed by a rotation, 
which can be written as (Ba§kal et al. 2005) 

5(e,A)i?(O,/i) = i?(0,Oi?M, (111) 

where 

cosh ^ = cosh fi cosh A + sinh fi sinh A cos 9, 

sin 6'[sinh A + tanh /i(cosh A — 1) cos 6] 



tan( 



sinh A cos 6 + tanh /u[l + (cosh A — 1) cos^ 6] 



2(sin6')[sinhAsinh/i + C_ cos^^] 
tan ^= ^^^^26) + 2 sinh A sinh cos ^ ' ^ ' 



with 



C± = (coshA±l)(cosh/i±l). (113) 
Indeed, Eq.f lllip can be written as 

R{co) = B{4>, -0 B{9, A) 5(0, /i), (114) 

three squeeze transformations lead to one rotation. 

We have done this calculation using the kinematics of Lorentz transforma- 
tions. On the other hand, it does not appear possible to do experiments using 
high-energy particles. However, it is gratifying to note that this experiment 
is possible in polarization optics. 

If the angle 9 is 90° we use the notation S'(A) for 5 (|, A). So 

S(A)=("*(^/2) -h(A/2)) 

Vsmh(A/2) cosh(A/2) / 

and the above calculations become simpler with 

Si\)B{0,fi)=B{(P,OR{co), (116) 
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where 



cosh ^ = cosh fi cosh A, 
sinh A 

tan ( 



tanh /i 
sinh A sinh /i 

tanu; = — — -. (H') 

cosh fi + cosh A 

6.3 Optical Activities 

For convenience, let us change the parameters 9 and n as 

e = 2az, fx = 2(3z, (118) 
and the R{0) and ^(/i) matrices as 

R{az) = C'^^'^'l -^i^M), (119) 
V sm(azj cos(azj / 

and the rotation angle increased as the beam propagates along the z direction. 
This version of optical activity is well known. 

In addition, we can consider the squeeze operation 

q( cosh(/3z) -sinh(/3z)\ , . 

^{ P^J- l^_sinh(/3z) cosh(/3^) )■ ^^^^> 

Here the squeeze parameter increases as the beam moves. The negative sign 
for /3 is for convenience. 

If this squeeze is followed by the rotation of Eq. flllQI) . the net effect is 



cosh(/32;) — sinh(/3z) \ / cos(a2;) — sin(az) 
— sinh(/3z) cosh(/3z) / v sin(a2;) cos(q;2;) 



121) 



where z is in a macroscopic scale, perhaps measured at the order of centime- 
ters. However, this is not an accurate description of the optical process. 

In fact it happens in a microscopic scale of z/N, and becomes accumulated 
into the macroscopic scale of z after N repetitions, where is a very large 
number. We are thus led to the transformation matrix of the form (Kim 
2010) 

M{a, (3, z) = [S{-(3z/N)R{az/N)f . (122) 
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In the limit of large A^, this quantity becomes 



1 -I3z/N\ ( 1 
-(3z/N 1 )\az/N 



-az/N 
1 



-, N 



(123) 



Since az/N and (5z/N are very small 

M{a,l3,z) = 
For large N, we can write this matrix as 



1 
1 







J N 



N 



M{a, (3, z) = exp (Hz) 



(124) 



(125) 



with 



H 



:i26) 



-{a + 13) 
(a-/3) 

We can compute this matrix using the procedure developed in Sec. [3l If 
a is greater than /3, H becomes 



\exp{—ri) 



(127) 



with 



a' = - /32 



exp (?]) 



Ia + (3 



I a — f3' 

and the M matrix of Eg. (11251) take the form 



(128) 



cos a z) 



e ' sm a z) 



e ^sm(a'z) cos(a'z) 



(129) 



If /3 is greater than a, the off-diagonal elements have the same sign. We 
can then write H as 



H = -(3' 







exp (r/) 



exp (—77) 



(130) 
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with 



and the M matrix of Eq. (11251) becomes 

(132) 



cosh(/3'z) —e^ sinh(/3'2;) 
-e"'' sinh(/3'z) cos\i{(5'z) 

If a = /3, the lower- left element of the H matrix has to vanish, and it 
becomes 

H-il (133) 



(134) 



,0 

and the M matrix takes the triangular form 

1 -2az 
1 



The optical material can be made to provide rotations of the polarization 
axis. It is much more interesting to see this additional effect of squeeze. 



7 Stokes Parameters and the Poincare Sphere 

In Sec. O we studied various aspects of the Jones vector formalism of the 
polarized beams, we have not dealt with the problem of whether the two 
beams are coherent with each other. In order to study this coherence problem 
we have to construct the four Stokes parameters. 
Let us write the Jones vector of Eq. ( IHOl) as 

U exp{z(fcz-u;t + 52)}y ' ^ ' 

where a and h are positive real numbers. In Sec. O we studied the effects of 
the squeeze B{rj), phase shift Z{5), and rotation R{0) on the Jones matrix. 
These matrices are given in Eq.f l83|) . Eq.f lHSj) and Eq. fl86|) respectively. 

These transformation matrices can be written as one expression as the 
two- by-two matrix of G of Eq. flBSl) . and its role in the Lorentz group and 
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its physical application to the Jones vectors were discussed in Sees. O and El 
respectively. While the Jones vector can deal with two independent beams, 
it does not address the issue of whether they are coherent with other. For 
this purpose, let us introduce the coherency matrix (Brosseau 1998, Saleh et 
ah 2007). 

C=(^'' ^''), (136) 

with 

< VrV-, >= ^ J%*{t + T)^j{t)dt, (137) 

where T is for a sufficiently long time interval, is much larger than r. Then, 
those four elements become 

Sn =< ^liJi >= a^, Si2 =< >= ab e~^''+'^\ 

S21 =< ij;iJi >= ab e-'^''-'^\ ^22 =< ^2>2 >= (138) 

The diagonal elements are the absolute values of ipi and ■?/'2 respectively. The 
off-diagonal elements could be smaller than the product of ipi and ip2, if the 
two beams are not completely coherent. Thus, the parameter a serves as the 
decoherence parameter. 

The a parameter specifies the degree of coherency. Unlike the ABCD 
matrix, this coherency matrix is not always real, and its determinant is not 
always one. If we restrict the trace of this matrix to be one by normalizing, 
this matrix becomes the density matrix (Feynman 1972). 

If we start with the Jones vector of the form of Eq.f ll35p . the coherency 
matrix becomes 

^-U&e-('^-^) b^ )■ ^^^^^ 



We are interested in the symmetry properties of this matrix. Since the trans- 
formation matrix applicable to the Jones vector is the two-by-two representa- 
tion of the Lorentz group, we are particularly interested in the transformation 
matrices applicable to this coherency matrix. 

The trace and the determinant of the above coherency matrix are 

det{C) = {abf (1 - , 

tr{C) = a^ + b^. (140) 
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Since e is always smaller than one, we can introduce an angle x defined as 

cosx = e-^ (141) 

and call it the "decoherrence angle." If x = 0, the decoherence is minimum, 
and it is maximum when x = 90°. We can then write the decoherency matrix 
of of Eq.(Il39D as 



C 



ab{cos x)e" 



,a6(cosx)e*^ 

The degree of polarization is defined as (Saleh et al. 2007) 



(142) 



P 



1 



4 det{C) 



1 



4(a6sin x)^) 



\ (tr(C))2 ^ (a2 + 62)2 • 
This degree is one if x = 0. It becomes 

a-' + b^' 



(143) 



(144) 



when X = 90°. We can without loss of generality assume that a is greater 
than b . If they are equal, the degree of polarization is zero. 



7.1 Stokes Parameters as Four- Vectors 

Under the influence of the G transformation given in Eq. f l65p , this coherency 
matrix is transformed as 



C' = G C 



or 

'-'11 




or 




Su 


S12 


S21 


S22 



a f3\fSu Si2\fa* 7* 
7 Sj\S2i S22J\/3* S* 



(145) 



Here, the G matrix is not Hermitian, and its Hermitian conjugate is not 
always its inverse. Thus it is not a similarity transformation, yet it pre- 
serves the determinant of G. When the G matrix in Eq.f l65|) consists of real 
elements, it becomes the transformation matrix applicable to the ABGD 
matrix. If it is constrained to be Hermitian, it becomes a rotation matrix 
without boosts. 
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While the coherency matrix is transformed as in Eq.f ll45p . its components 
transforms in the same manner as in Eg. fl67l) 



(S'u\ 




/ aa* 


a(3* 


(3a* 


•^12 




aY 


aS* 


h* 


CI 

'-'21 




* 


7/3* 


6a* 






\77* 


75* 


57* 



/3/3*\ 
136* 
5/3* 



\ 


(Sii\ 




S12 




S21 


/ 


\S22/ 



(146) 



Particular combinations of the coherency matrix components of Eq.f ll36p 
are crucial for the quantum picture of polarization (Falkoff et al. 1951): 



So 



Si 



Su + "S*: 



22 



V2 ^ 

5*12 + 5*21 



V2 



5, 



^2 



Su — S; 



22 



V2 ^ 

5*12 — 5*21 



V2i 



(147) 



which can also be expressed as the sum of an identity matrix and the Pauli 
spin matrices a* (Fano 1954) 



C = -{Sol + Sia^ + S2a^ + S^a^) 



(148) 



where the coefficients are known as the Stokes parameters in the litera- 
ture (Shurcliff 1962). Furthermore, expressing in such a compact form also 
serves to treat the Jones and Mueller calculi in the framework of pure oper- 
ational Pauli algebraic approach (Tudor 2010). 

The the four-by-four matrix which transforms {Su, S22, S12, S21) to (5*0, 5*3, Si, S2) 

is 



(149) 



This matrix enables us to construct the transformation matrix applicable to 
the Stokes parameters, widely known as the Mueller matrix (Soleillet 1929, 
Mueller 1943, Brosseau 1998). 

The Mueller matrix applicable to the Stokes parameters takes the same 
form as the Lorentz transformation matrix applicable to the space-time four- 
vector of {t, z, X, y) given in Sec. [5l Therefore, the Mueller matrix is a four- 
by-four representation of the Lorentz group. 



fSo\ 






1 







fSiA 


S3 


1 


1 


-1 








S22 


Si 


"71 








1 


1 


S12 


\S2) 




\o 





—i 




\S2j 
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It is gratifying to note that the four-by-four Mueller matrices share the 
same symmetry properties as those of the two-by-two Jones matrices appli- 
cable to the Jones vectors. Thanks to the squeeze matrix -B(/i) of Eq. fj83|) . we 
are able to extend the symmetry of those two-by-two matrices from SU{2) 
to SL{2,c) (Han et al. 1997, Devlaminck et al. 2008, Redkov 2011). 

We should note here that the decoherence angle is a Lorentz-invariant 
quantity. It cannot be changed by Mueller transformations. It may be 
possible to construct a four-by-four matrix which will change this param- 
eter (Ortega-Qujiano et al. 2011), but this matrix cannot belong to the 
Lorentz group. 

As for the two-by-two matrix formalism, it is an iterating proposition to 
formulate the problem using quarternions (Dlugunovich et al. 2009, Tudor 
2010). Quarternions represent a four-dimensional rotation group and do 
more. Thus, interesting results may be obtained from this line of approach. 



7.2 Winger's Little Group for Internal Space- time Sym- 
metries 

It is more interesting to study the problem using the two-by-two representa- 
tion of the coherency matrix because their elements are directly measurable 
quantities. For this purpose, let us recall the Lorentz transformation of the 
four- vector of a free particle (po,P3,Pi,P2) which is the same as that of the 
{t, z, X, y) four-vector, and we can write its G transformation as 

P' = G P = ( P'-''P'y 



7 ^ )\Px + iPy Po-Pz J \P* S* 



(150) 



The two- by-two matrix for the four- momentum P is given in Eq. fl71l) . 

We can consider transformations which will leave the four-momentum 
invariant. In other words, we can write Wigner's little group as the subset 
of the G matrix which satisfies 

fPo+Pz Px-iPy\^fcy P\fPo+Pz Px-iPy\fa* 7*\ 

\Px + iPy Po-Pz) V7 SjKPx + iPy Po - Pz J \ 1^* S* J ' ^ 

Using the rotation matrix Z{6) of Eq. (!85|) which leads to a phase shift, we 
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can bring this formula to the form 

fpo+Pz Px \ ^ f a (3\ fpo+Pz Px W a 7 
V Px Po-Pz) V 7 ^ / V Px Po-Pz) \/3 S 

This can then be transformed to a diagonal form 



po + \pi+pi 



(152) 



(153) 



po-^pI+pI 

with the rotation matrix R{C)) where 

tan^ = — . (154) 

Pz 

With the boost squeeze matrix B{ri) given in Eq.( l83l) where 



= (155) 

\ Po-Pz 



we can transform the diagonal matrix of Eq. f ll53p to another diagonal matrix 

(156) 



m 

m 

where m = ^Jpl — p1~ p^- The squeeze matrix B{ri) corresponds the Lorentz 
boost given in Eq.f l72|) . 

The G transformation matrix which will leave this four-momentum matrix 
is has to be Hermitian, and is a rotation matrix R{6). This defines Wigner's 
internal space-time symmetry of a massive particle. 

This matrix remains invariant under the G transformation if the G matrix 
is Hermitian. It is Hermitian only for rotations. This is thus consistent with 
Wigner's 0(3)-like little group for massive particles, as discussed in Sec. [51 

For a massless particle, we can choose the system where pq = p^ = u and 
Px = Py = ^- Then the P matrix becomes 

and its determinant is zero, saying p^—pl = 0. It is not difficult to construct 
the G matrix whose G transformation will leave the above P matrix invariant. 
It takes triangular form 

1 7' 
1 



(158) 
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From this, it is not difficult to construct its four-by-four counterpart given 
in Eq.dlH]). 

Let us go back to the symmetry of the coherency matrix, and to the 
matrix C of Eg. (11391) . If we make a G transformation with Z{5), the D 
matrix becomes 

a? ah cosx~ 
.a6 cosx 

The two-by-two matrix of Z{5) is given in Eq.( 

If we make another G transformation with R{—0) where 



(159) 



2ab cosx .^„„^ 

tan6' = — — , (160) 

a'' — 0^ 



the coherency matrix becomes 

J' 

s — r 



C={'t'' ° ), (1611 



with 



' 2 

r = {a? + h^f + A{ahY sin^ X- (162) 

If ^ = 0, the system is totally coherent, and the coherency matrix becomes 

-a2 + 62 0- 







(163) 



This matrix is like the four- momentum matrix of Eq. fl7T]) for massless parti- 
cles. 

If X is nonzero, we can G-transform the D matrix of with 

B{-V)={\ ^'J, (164) 



with 

Is + r 



s — r 



the D matrix becomes 
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This coherency matrix is invariant under G transformations if the G matrix 
consists only of rotations and thus is Hermitian. This aspect is consistent 
with Wigner's 0(3)-hke httle group for massive particles. 

In the case of the four-momentum matrix, its determinant is m? and is 
Lorentz-invariant. For the coherency matrix, the determinant is (a&)^ sin^ x- 
This means that the coherency angle x cannot be changed by Lorentz trans- 
formations, as in the case of mass in special relativity. 



8 Symmetries of the Poincare Sphere 

The Poincare sphere has a long history, and its spherical symmetry is well 
known (Born et al. 1980). The rotational symmetry of the Poincare sphere is 
readily included in the Lorentz symmetry. We shall first review the rotational 
symmetry, and study the effect of Lorentz boosts. 

Let us write the coherency matrix of Eq.f ll39p as a four-component vector 







( ^ 


\ 


rz 




r(cos 


0) 






r (sin 6') 


cos 5 


Vy) 




\ r(sin^) 


sin 5 / 



/ (a2 + 62)/2 
(a2-62)/2 
a6(cos 5) cos x 
V a6(sin 5) cosx / 



(166) 



This four-vector is defined by two spheres: the outer sphere with radius s, 
where 

.^<^. (167, 

and the inner sphere defined by the three-component vector (r^, r^^, r^), with 
its radius 



= {a? — 6^)^ + 4(a6)2 cos^ x 
which is the radius of the Poincare sphere. Its z component is 



(168) 



(169) 



which is independent of the decoherency angle x- Here we assume the am- 
plitude a to be greater than h. 
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The radius of the Poincare sphere r depends on x, and takes its maximum 
value s when x = 0. The radius shrinks to its minimum value when % 
goes to it largest value. Figure H] illustrates the circles and their radii. The 
ratio of to r is 

cosg = ^ = ^ (170) 
If we apply the rotation 



/I 













1 














COS 5 


sin 5 


Vo 





— sin 5 


cos 5 j 



(171) 



and then 



/I 













COS 6* 


sin Q 








— sin Q 


cos Q 





vo 








1/ 



(172) 



with cos 6' = r^/r given in Eg. (11 70 p . this four- vector can be brought into the 
form 

(s, r, 0, 0). (173) 

Within the framework of the traditional three-dimensional geometry of 
the Poincare sphere, it is possible to transform the four vector of Eq.f ll66p 
to the four vector of Eg. (11 73 p while the value of s is left unchanged. On 
the other hand, the Lorentz symmetry allows a transformation on this four- 
vector by B{—ri) with tanhi] = r/s so that the four- vector becomes (Ba§kal 
et al. 2006) 





^2 \ 


/ 


cosh rj 


— sinh rj 
















— sinh 1] 


cosh 7] 








r 















1 








V 


/ 


V 











1/ 


VO/ 



(174) 



This means that the radius of the Poincare sphere can become zero while the 
outer radius takes its minimum value ab sinx- 

Indeed, the angle x determines the minimum radius of the outer sphere. 
This radius takes the maximum value of ab when a becomes infinity. The 
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larger radius becomes zero when a = corresponding to completely coherent 
beams, and it does not correspond to the real world. 

In order to resolve this problem, we note that the large radius and smaller 
radius are the same and its value is (a^ + 6^) /2, when x = 0. Thus, if we 
make the Lorentz boost of Eq. f ll74p , the result is 





/ 


s 











V 



cosh rj — sinh rj 

— sinh 1] cosh rj 

1 





\ 






s 












(175) 



The Lorentz group, including its rotation subgroup, changes all the pa- 
rameters for the coherency matrix. However, it cannot change the decoher- 
ence angle x- What significance does this carry from the symmetry point of 
view and from the physical point of view? 



Outer Sphere 




Figure 4: Poincare sphere and its outer sphere. The radius of the Poincare 
sphere depends on the decoherence angle x, but its z component does not. 



8.1 0(3,2) Symmetry of the Poincare Sphere 

It is clear from the previous subsection that the decoherence angle x is invari- 
ant under Lorentz transformations. We are now looking for another symme- 
try which will change this variable. For this purpose we write the coherency 
matrix of Eq. f ll42p as 

_ / ab{cosx)e-'^\ . „ . 

[ab{cosx)e'' J ^ ' 
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and introduce another matrix where cosx is replaced by sinx, which can 
take the form 

a6(sinx)e~ 
, a6(sin x)e*'' 6^ 

For the first coherency matrix of Eq. fll76p , we know how to transform it 
into the four- vector (s, ri, 0, 0), with 



(177) 



n = I ^J{a^ - b^y - A{aby {cos xY- 



(178) 



Then, the second matrix can be brought to the four- vector (m, r2, 0, 0), 
with 

u = ^-^, r2 = ^V'(a2-62)2_4(a5)2(sin^)2. (179) 
These two expressions lead us to consider the following two four vectors. 



(180) 



In view of the relation cos^ x + sin^ x = 1, we are led to the five dimensional 
vector space with (s, u, r^, r^, Ty), which can start with 



/ a6(cos x) \ 




/ a6(sin x) \ 


















V ) 




I / 



icosx, afosinx, 0, 0, 0) 



;i8i) 



Now we can change the value of the decoherence parameter a by changing 
the angle x? but we can change this variable by introducing a rotation matrix 
applicable to the two-dimensional vector space of s and t. When all other 
components vanish we can write the rotation matrix 



ab{cos x) 
a6(sin x) 



cosx 
sinx 



sin x\ ( (^b 



cosx 7 V 



182) 



Therefore the five-by-five rotation matrix will be of the form 

/ cos X — sin X \ 

sinx cosx 

10 

1 

V 1/ 



(183) 
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which is apphcable to the five-component vector (s, t, r^, r^, Vy), with the 
two four-dimensional subspaces, corresponding to the coherency matrices 



Csix) 



a? abe *^(cosx) 

6ae*''(cosx) b"^ 



Here both — and — remain invariant under Lorentz transformations 
within their four-dimensional subspaces. Their deteminants are 

det{C s) = (abY sin^ X, det{C u) = {abf cos^ x, (185) 

resuting in 

det (Cs) + det (Cs) = (abf, (186) 
which is independent of the angle x- 

8.2 Entropy Probem 

Another important way to measure the lack of coherence is to calculate the 
entropy of the system. The coherency matrices defined in this section become 
density matrices if their traces are normalized to be one. Thus from Eq. fll77p . 
we can write 



Psix) 



abe (cosx) 



_|_ ^2 \^ la e*'^(cos x) b' 



_ 1 f a^ abe ^'^(sinx)\ 

^-^^^-^^^Uae^^(sinx) b' 



These matrices can be diagonalized to 

Psix) = ^ ( 



^fl + fs 

1-/. 



where 
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Then, their entropies become 



The entropy Sg becomes zero when x = 0. It becomes 

when X = 90°. The entropy 5*^ is a monotonically increasing function of the 
angle x starting from zero to the above maximum value, which becomes In 2 
when a = b. 

The entropy Sg of the first space is monotonically increasing function of x, 
while that of the second space Su is a decreasing function. Thus, an increase 
in entropy in the first space leads to a decrease in the second space. Then 
we can ask whether the sum of these two entropies becomes independent of 
X, leading to an entropy conservation of the total system. The answer is No. 
However, this does not cause problems for us, because the second space is 
not necessarily a physical space. It could be meaningless to use the same 
definition of entropy for this space. On the other hand, as we noted before, 
we can define the conservation of entropy in terms of the the sum of the 
determinants of the coherency matrices given in Eq.f ll84p . Furthermore, this 
determinant condition does not require that the amplitudes of the two beams 
to be the same. 

8.3 Feynman's Rest of the Universe 

What is the meaning of this second space? In his book on statistical mechan- 
ics, Feynman makes the following statement about the density matrix (Feyn- 
man 1972): When we solve a quantum-mechanical problem, what we really do 
is divide the universe into two parts - the system in which we are interested 
and the rest of the universe. We then usually act as if the system in which we 
are interested comprised the entire universe. To motivate the use of density 
matrices, let us see what happens when we include the part of the universe 
outside the system. 
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Feynman did not specify whether the rest of the universe is observable 
or not. In either case, it is an interesting exercise to construct a model of 
the rest of the universe behaving like a physical world. With this point in 
mind, one of us studied two coupled harmonic oscillators in which one of the 
oscillators correspond to the physical world and the other to the rest of the 
universe (Han et al. 1999). In this example, the rest of the universe is the 
same as the world in which we do physics. In thermal field theory (Umezawa 
ct al. 1982), even though based on the same mathematics as that of the 
coupled oscillators, the rest of the universe is not physically identified, except 
that it causes thermal excitations of the oscillators in the physical world. 

The concept of decoherence occupies one of the central places in the cur- 
rent development of physics. In (Feynman et al. 1963) when the system 
couples to finite temperature baths the result is an environmentally induced 
decoherence. Its effects can be determined in tunneling processes (Caldeira 
et al. 1983) and in two-state systems that are coupled to dissipative en- 
vironments (Leggett et al. 1987). The decoherence of the electromagnetic 
field coupling resonantly to a two-level system (Anastopoulos et al. 2000) 
and applications of two-level decoherence to qubit systems (Shiokawa et al. 
2004) are also investigated in the literature. 

The pattern for the two-optical beams arising from phase-randomizing 
process (McAhster et al. 1997), is in the same structure of the two-by-two 
matrix discussed in this paper. As for the decoherence in the rest of the 
universe introduced in this work, the system becomes more coherent as the 
time- variable increases. Although this "recohcrcncc" process was considered 
earlier in the literature (Anglin et al. 1996), it is premature to expect a 
two-state system to gain coherence in the real world. It is thus very safe to 
say that the second Minkowskian space introduced in this paper remains in 
Feynman's rest of the universe. 

However, this does not prevent us from constructing a physical system 
analogous to the decoherent system coupled to a recoherent system. 

Concluding Remarks 

In this paper, we have organized ray and polarization optics using the lan- 
guage of the Lorentz group. The Lorentz group has two-by-two and four- by- 
four representations. Both are useful in optics, and they allow us to gain a 
unified view of various aspects of optics. 



42 



In addition, it was noted that the mathematics apphcable to ray and po- 
larization optics is directly applicable to the internal space-time symmetries 
of elementary particles. Optical systems are favorable in the sense that each 
mathematical operation has a counterpart that can be performed in optics 
laboratories. 

The Lorentz group is also the basic mathematical language for Einstein's 
special relativity. Currently, this group serves useful purposes in many other 
branches of physics, including optical sciences. In recent years, the Lorentz 
group served as the underlying language for squeezed states of light. It 
was Dirac who first observed that the Lorentz boost is a squeeze transforma- 
tion (Dirac 1949) and constructed representations of the Lorentz group using 
coupled harmonic oscillators (Dirac 1963). Indeed, Dirac's oscillator repre- 
sentation forms the theoretical foundations of squeezed states of light (Yuen 
1976, Yurke et al. 1986, Kim et al. 1991). This aspect of the Lorentz group 
is by now well known in the optics community, and the Lorentz group is one 
of the theoretical tools in quantum optics. 

The squeezed state is not the only branch of optics requiring the Lorentz 
group. It can well be applied to Fourier optics (Bacry et al. 1981) while its 
geometry has proven to be useful in designing three dimensional non-imaging 
concentrators (Gutierrez et al. 1996). Para-axial wave optics (Sudarshan et 
al.l983, Makunda et al. 1983) and Wavelets (Han et al. 1995) are also 
known to be representations of this group. It is also the underlying language 
for reflections and refractions (Pellat-Finet et al. 1992). 
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